Astronomy & Astrophysics manuscript no. cepatherm'27a 


©ESO 2013 


January 12, 2013 





Distribution and excitation of thermal methanol in 6.7 GHz maser 

bearing star-forming regions 

I. The nearby source Cepheus A 

Karl J.E. Torstensson 12 , Floris F.S. van der Tak 3 ' 4 , Huib Jan van Langevelde 21 , Lars E. Kristensen 1 , and Wouter H.T. 

Vlemmings 5 

1 Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands 
e-mail: [kalle ; kristensen] @strw . leidenuniv . nl 

2 Joint Institute for VLBI in Europe, PO Box 2, 7990 AA Dwingeloo, The Netherlands 
e-mail: langevelde@jive.nl 

3 SRON Netherlands Institute for Space Research 

C^i [ Landleven 12, 9747 AD Groningen, The Netherlands e-mail: vdtak@sron.nl 

D ■ 4 Kapteyn Astronomical Institute 

■ University of Groningen, The Netherlands 

5 Argelander-Institut fur Astronomie, University of Bonn 
Auf dem Hugel 71, D-53121 Bonn, Germany e-mail: wouter@astro.uni-bonn.de 



o 

(N 



oo 



in 



On 



in 



oo 



- 1—1 

x 



Received September 15, 1996; accepted March 16, 1997 

ABSTRACT 



Context. Candidate high mass star forming regions can be identified through the occurrence of 6.7 GHz methanol masers. In these 
sources the methanol abundance of the gas must be enhanced, as the masers require a considerable methanol path length. The place 
and time of origin of this enhancement is not well known. Similarly, it is debated in which of the physical components of the high 
Q ' mass star forming region the masers are located. 

Aims. The aim of this study is to investigate the distribution and excitation of the methanol gas around Cep A and to describe 
the physical conditions of the region. In addition the large scale abundance distribution is determined in order to understand the 
C$ ' morphology and kinematics of star forming regions in which methanol masers occur. 

Methods. The spatial distribution of the methanol is studied by mapping line emission, as well as the column density and excitation 
temperature, which are estimated using rotation diagrams. For a limited number of positions the parameters are checked with non-LTE 
models. Furthermore, the distribution of the methanol abundance is derived in comparison with archival dust continuum maps. 
Results. Methanol is detected over a 0.3 x 0.15 pc area centred on the Cep A HW2 source, showing an outflow signature. Most of the 
gas can be characterized by a moderately warm rotation temperature (30-60 K). At the central position two velocity components are 



detected with different excitation characteristics, the first related to the large-scale outflow. The second component, uniquely detected 



at the central location, is probably associated with the maser emission on much smaller scales of 2". Detailed analysis reveals that the 
highest densities and temperatures occur for these inner components. In the inner region the dust and gas are shown to have different 
| physical parameters. 

Conclusions. Abundances of methanol in the range 10~ 9 - 10~ 7 are inferred, with the abundance peaking at the maser position. The 
geometry of the large-scale methanol is in accordance with previous determinations of the Cep A geometry, in particular those from 
methanol masers. The dynamical and chemical time-scales are consistent with the methanol originating from a single driving source 
associated with the HW2 object and the masers in its equatorial region. 
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1. Introduction Methanol maser sources are exciting targets for the de- 
tailed study of high-mass star-formation. Continuum emission 

While massive stars (>1OM ) have an enormous impact on the of warm dust at sub . millime t r e wavelengths has been detected 

evolution of the Galaxy, during their life and death their birth at well ove r 95% of the observed 6.7 GHz methanol maser sites 

process is still shrouded in mystery (e.g. | Kurtz || 2005|) . To a large dHill et al]l2005l) . while 25% of the same masers ar e associated 

extent this is due to the environment in which the high-mass stars with a detectable u it ra -compact (UC) Hn region dWalsh et al.l 

form. Most high mass stars are observed at large distances (com- 1998) This seems tQ indicate ^ ^ mafjers probe a range of 

pared to low mass stars), evolve rapidly, and form deeply embed- early phases of massive star formation and that the met hanol 

ded in their molecular clouds. Additionally, they typically form maser emission disa pp ears when the UC Hn region is created 

in clusters. Thus, very high resolution data are needed to resolve ( WaM et al 199g) The maser sources require high methanol 

them and the disruptive influence they have on their direct envi- abundance for meir excitatio n, pointing to a re cent release of 

ronments (e.g.| Martin-Pintado et al. || 2005 | ) . methanol into the gas phase dvan der Tak et al.ll2000h . We can 

thus expect to observe the thermal methanol lines in these 
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Fig. 1. Sample spectra from the 3 positions identified in Sec. 13.21 Top spectrum at the "Centre" position, coincident with the HW2 
object. Bottom Left at the position of "NE Outflow" and Bottom Right "Envelope". The small symbols at the top indicate the 
frequencies of methanol transitions in this band. 



sources, yielding diagnostics for the large scale kinematics and 
excitation. Besides being signposts, masers are excellent trac- 
ers of the geometry and small scale dynamics of these regions. 
Much effort has focused on studying the kinematics and sev- 
eral claims of circumstellar di sks, expanding spherical shells and 
jets have been presented (e.g. [ Pestalozzi et al. 2004; Minier et~aD 
l2002tlBartkiewicz et al.l2005ll2009t) . The evidence suggests that 
masers trace an evolutionary sequence or possib ly a mass rang e 
in the formation process, or even both (e.g. iBreen et al.l l2010). 
However, it is not clear exactly with which of the physical struc- 
ture^) of the complex surroundings of young massive stars they 
are associated. By trying to make a link between the maser dis- 
tribution and the large scale methanol kinematics and excitation, 
one can attempt to reach a better understanding of the physical 
and chemical state of the gas of the regions in which masers oc- 
cur. 

In order to carry out such a study, we have observed thermal 
methanol gas towards a sample of 15 high-mass star-forming 
regions for which we have very long baseline interferometry 
(VLBI) data of the 6.7 GHz methanol maser. The observations 
yield maps of the distribution and physical condition of the 
methanol gas in relation to the methanol maser. For one, the 
observations of the methanol allow one to probe a reasonably 
large range of temperature and density, due to its closely spaced 
AA' = (A7 =1) transitions with a large range i n excitation en- 
ergies and critical densities dLeurini et al.l 12004). Moreover we 



can measure the methanol abundance and velocity field. Clearly, 
the metha nol abundance need s to be significant in these regions 
(> 10~ 7 . ISobolev et alJll997l) : evaporation of CH 3 OH from the 
dust grains in these regions, at temperatures ~100 K is thought 
to be responsible for this. This specific enhancement is supposed 
to be s o short lived that it can be used as a chemical clock of the 
region (Ivan der Tak et al.ll2000h . The enhancement can be used 
to determine the age of these regions and thus the evolutionary 
stage of the star formation process. 

In this paper, we describe the observations and the procedure 
we have used to analyse the data. Because of its proximity and 
because a wealth of data exists on this source, we will focus on 
Cepheus A East as a first example. The other sources will be 
presented in a second paper. 

At a distance of 700 pc dMoscadelli et alJl2009h . Cepheus 
A East is one of the closest high-mass sta r forming regions . 
The region has a luminosity of 2.3 x 10 4 L Q (lEvans et alJ ll981) 
(scaled to the distance of 700 pc), out of w hich half is though t 
to come from the well-studied object HW2 (Ga rav et al.l ll996). 
which is also the site of the 6.7 GHz methanol maser emission. 
The object shows b oth large scale molecular outflows over ~1' 
(Gomez et al. 1999) a nd signs of collapse over a similar extent 
(Botfinelli & Williams! 120041: ISun & Gadl2009l) . At the core of 
HW2, a thermal jet extends over ~2" and proper motion out- 
flow velocities of ~500 kms' 1 have been measured in the ra- 
dio dHughes & Wouterlootl[l9"84t ICuriel et al.ll2006h . Although 
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Table 1. CH3OH line data for the observed transitions, adopted 
from the CDMS (C ologne Database of Molecular Spectroscopy, 
iMtiller et ail (12005b '). All lines are from the J = 7 -> 6 band, 
and throughout the paper a notation like "-6 E" refers to the 
J = 7_e, — > 6_6 E transition. Blended lines are indicated by a *. 



Frequency 


*v g 


E„ 


Transition 


MHz 


D 2 


K 


k type 


337969.414 


5.55 


390.1 


-1 A v,= l 


338124 502 


5.65 


78.1 


E 


338344.628 


5.55 


70.6 


-1 E 


338404 580 


1.49 


243.8 


+6 E 


338408.681 


5.66 


65.0 


A 


338430 933 


1.50 


253.9 


-6E 


338442.344* 


1.49 


258.7 


+6 A 


338442.344* 


1.49 


258.7 


-6 A 


338456.499 


2.76 


189.0 


-5 E 


338475 290 


2.76 


201.1 


+5 E 


338486.337* 


2.77 


202.9 


+5 A 


338486 337* 


2.77 


202.9 


-5 A 


338504 099 


3.80 


152.9 


-4 E 


338512.627* 


3.81 


145.3 


-4 A 


338512.639* 


3.81 


145.3 


+4 A 


338512.856* 


5.23 


102.7 


-2 A 


338530.249 


3.82 


161.0 


+4E 


338540.795* 


4.60 


114.8 


+3 A 


338543.204* 


4.60 


114.8 


-3 A 


338559.928 


4.64 


127.7 


-3E 


338583.195 


4.62 


112.7 


+3E 


338614.999 


5.68 


86.1 


+ 1 E 


338639.939 


5.23 


102.7 


+2 A 


338721.630 


5.14 


87.3 


+2E 


338722.940 


5.20 


90.9 


-2E 



a superposition of gas comp onents at different positions can- 
not b e completely ruled out dBrogan et al.l I2007T IComito et alJ 
120071) . the evidence is building that there is a disk a t the centre 
of the jet, observable in both molecular gas and dust dPatel et all 
12005k Torr elles et al.l2007tlJimenez-Serra et al.l2007T) . as argued 
in detail in Ijimenez-Serra et all (120091) . VLBI observations of 
the 6.7 GHz methanol maser suggest that the methanol masers 
arise in a ring-li ke structure, extending over ~2", straddling the 
waist of HW2 dSueivama et al.l I200I IVlemmings et all 12010k 
Torst ensson et al.ll201 ll) . 

In Sec. [2] we describe the observations and calibration meth- 
ods. The results are presented in Sec. [3] in the form of a num- 
ber of sample spectra and maps of integrated line strength and 
velocity. In order to study the large scale methanol excitation 
we construct rotation diagrams and produce maps of the derived 
physical quantities in Sec. [4] We present non-LTE calculations 
for a number of positions to evaluate the limitations of the di- 
agram analysis in Sec. 14.41 Subsequently we discuss estimates 
of the hydrogen column density before we present the methanol 
abundance in Sec. 14.61 The relation of our findings to the Cep 
A star forming process and in particular the occurrence of the 
methanol maser are discussed in Sec. [5] 

2. Observations and data reduction 

The JCMlQ observations of Cep A East (HW2) (22 /! 56 m 17.9 ! 
+62°01'49" (J2000)) were performed on June 2 and June 17, 

1 The James Clerk Maxwell Telescope is operated by the Joint 
Astronomy Centre on behalf of the Science and Technology Facilities 
Council of the United Kingdom, the Netherlands Organisation for 
Scientific Research, and the National Research Council of Canada. 



2007. We used the array receiver HARP, which has 16 receptors 
with a spacing of ~30" in a 4x4 grid, resulting in a footprint of 
~2'. The receptors are single side-band receivers with a sideband 
rejection of >10 dB. We used the observing mode HARP5, a 
type of beam-switching (jiggle-chop) mode where several short 
scans at different positions of the target source (10 s per jiggle 
position) are observed before switching to the off position. As a 
result, several on-source observations share the same off-source 
observation which minimises the overhead while still providing 
good baselines. The observing mode results in a map of 20 x 20 
pixels covering 2'x2', thus each pixel is 6"x6". This ensures 
proper Nyquist sampling of the 14" telescope beam at 338 GHz. 

The ACSIS correlator backend was set up with a 1 GHz 
(880 kms -1 ) bandwidth and 2048 channels centred on the 
methanol 7o — > 60 A + line at 338.41 GHz. The frequency set 
up allows us to cover 25 CH3OH lines (Table [TJ with a veloc- 
ity resolution of ~0.43 kms" 1 (488 kHz). In total we have an 
effective (on source) integration time of ~5 min per pixel with 
a typical system temperature (r sys ) of ~280 K. For our analysis 
we have adopted a main beam efficiency (^mb) of 0.6. Regular 
pointings were done of the normal calibrators and we estimate 
to have an absolute pointing accuracy of «sl". 

We used the Starlink package (Gaia and Splat) for initial in- 
spection of the data. It was found that receptor R06 suffered from 
bad baselines (large ripples is some scans). For that receptor only 
these poor scans were removed, after which the scans were com- 
bined. The data were then converted to GILDAS/CLASS format 
and the remaining data reduction and analysis was performed in 
CLASS. A linear baseline was fitted to several spectral regions 
without any line emission and subtracted from the spectra. Next, 
to improve the signal to noise, the data were smoothed in the 
spectral domain resulting in a velocity resolution of 0.87 km s _1 
and an rms of ~ 65 mK. 

The analysis was performed on a pixel by pixel basis, by first 
fitting a Gaussian to the strongest unblended line (7_i — » 6_i E) 
in the spectra. The measured velocity and width of this line was 
then used to place windows around all other features in the spec- 
tra, for which the emission was integrated. Due to some overlap 
between different lines, only lines with a > 5<x (~ 0.6 Kkms -1 ) 
detection are included in the subsequent analysis. In the centre 
area, at the position of HW2, a second velocity component can 
be distinguished. Where the two velocity components could be 
separated, manual analysis by fitting individual Gaussians to the 
lines was performed. 

3. Results 

3.1. Methanol lines 

In the following we refer to each pixel in the map with its re- 
spective offset coordinates (Aa, Ad), measured in arc seconds, 
with respect to the centre of the map (J2000 22 /, 56'" 17.88* 
+62°01'49.2"). Fig. [TJ shows the observed spectrum at 3 posi- 
tions: (+3",3") "Centre", (+27",+9") "NE outflow", (+21",-3") 
"Envelope", where the designations can be understood from the 
forthcoming discussion in Sec. 13.21 Focusing at the centre spec- 
trum, all but four of the detected lines are due to methanol (indi- 
cated by ticks at the top); the others are identified as SO2, 34 S02 
and H 2 CS. The CH3OH K=+l E line is a blend with another 
SO2 line, also indicated in the figure. 

In the centre spectrum, two gas components are observed: 
the main component with Vi sr = -10.5 km s~' and a secondary 
red shifted at — 5.2kms _1 , as can be seen in the profile of the 
-IE transition (Fig. |2j. We note that the red-shifted velocity 
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Fig. 2. Velocity profiles of the - IE line at 3 positions. Although 
the profiles are Gaussian in the outflow, there is a clear shift in 
the velocity. At the "Centre" position there are clearly two veloc- 
ity components present. The dotted line is taken to be the system 
velocity, corresponding to the velocity of the dominant compo- 
nent at the "Centre" position. 



profile is overlapping with the range spanned by the methanol 
masers, which a re observed with veloc ities between vi sr = -4.6 
and -1.5kms _1 dTorstensson et al.ll201 lb . The red-shifted com- 
ponent is most readily identified in the higher-K transitions. 
While the first component is brightest in the low-K transitions 
and not detected in the higher-K lines, the — 5.2 kms compo- 
nent is brighter in the higher-K lines and less bright in the lower- 
K lines. For the K=3 transitions, the two components are of sim- 
ilar strength. There is a weak detection of the methanol isotopo- 
logue 13 CH 3 OH J = 13 -» 12, Ak = +1 in this component, 
indicating optical depth effects may be important for this. Also 
careful analysis, including averaging over a few pixels, shows 
that the torsionally excited -1 A v t =l line at 337.969 GHz is 
detected in the "Centre" region in the -5.2 kms 1 velocity com- 
ponent. 

Otherwise the observed line profile appears Gaussian in 
all components and does not show evidence of blue/redshifted 
wings. However, a modest velocity shift of the line is observed 
across the source with lower velocities to the NE and higher ve- 
locities to the SW (Fig.©. 

In order to start a consistent analysis of all methanol line 
emission in our data cube, we have extracted the integrated line 
flux from each line as described above (Sec. [2). To illustrate our 
procedure, we list the integrated line emission of the spectral fea- 
tures for the main velocity component at the "Centre" position 
from Fig.[T]in Table|2] For lines that are blended we have chosen 
a simple strategy. If the strengths (S g fig) of the lines are similar 
and they have comparable upper energy levels, we have split the 



Table 2. Line parameters for the main velocity component mea- 
sured at the position of HW2 (+3", +3") 



Line 


Flux 


Width 


K, type^ 


Kkms 1 


kms" 1 


0E 


2.16 


4.21 


-1 E 


3.30 


4.79 


OA 


3.60 


4.79 


-6E 


0.03 


1.74 


+6 A 


0.29 


3.80 


-6 A 


0.29 


3.80 


-5E 


0.61 


4.81 


+5E 


0.31 


3.15 


+5 A 


0.36 


4.18 


-5 A^ 


0.36 


4.18 


-4E~ 


1.05 


4.79 


-2 A* 


2.00 


4.79 


+4E~ 


1.12 


4.79 


+3 A 


1.37 


4.79 


-3 A^ 


1.37 


4.79 


-3E~ 


0.93 


4.07 


+3E 


0.94 


4.97 


+2 A 


1.42 


5.38 


+2E 


1.55 


4.38 


-2E" 


1.55 


4.38 



" See Table 1 for transition details 
b Blend with +5 A 
Blend with +4 A 
d Blend with +3 A 
' Blend with +2 E 



measured flux between the two lines. If however one line has a 
significantly higher line strength and lower upper energy level 
we attribute all the emission to the strong line. In the case of the 
blend between the K=-2 A and K=±4 A, where the lines have 
comparable line strengths and energy levels we assume the total 
measured flux as an upper limit to all the lines. 

3.2. Spatial distribution of methanol 

Fig.[3]shows the spatial distribution of the flux, the central veloc- 
ity, and the width of the methanol 7_i — > 6_i E line, which is the 
brightest unblended feature in the spectra. The methanol maser 
emission arises just to the NE of the centre of the map (+0.7", 
+0.4") and is indicated with a cross in all the figures. The white 
area in the maps represents regions where the line flux is less 
than our 5cr detection limit of 0.6 K km s _1 . 

The line emission is constrained to an area, measured from 
the points furthest separated in the intensity map (Fig.|3^), with 
an extent of -85" by 43" (-0.29 pc by -0.15 pc), and is ori- 
ented in the NE-SW direction. The velocity map (Fig.[3j3) shows 
that the methanol gas is entrained in an outflow with the blue- 
shifted emission to the NE and the red-shifted emission to the 
SW. Measured between the highest and lowest velocities in the 
map, the outflow has a modest velocity gradient of ~3.1 kms" 1 
over -0.20 pc. The integrated line flux peaks in the centre of the 
map at - 5.2 Kkms" 1 , with a bright extension to the SW over 
~13"(~0.044 pc), which is associated with the red-shifted emis- 
sion. Both the velocity field and line width maps (Figs. [3t> and 
c) show a maximum in the centre of the map to the N of HW2. 
This is due to blending of the two velocity components that we 
are unable to decompose. A second maximum in the integrated 
line flux is observed ~28"(~0.095 pc) to the NE of the centre, 
which is associated with the blue shifted part of the emission and 



Karl J.E. Torstensson et al.: Thermal methanol in Cep A 5 



x 
E 

2 -10 



3 
x 

E 




50 40 30 20 10 -10 -20 -30 -40 -50 
offset from HW2 (E-W) ["] 



50 
40 

i 20 

1 o 

E 

g -10 

1 -20 
o 

-30 

-40 
-50 




* * 



50 40 30 20 10 -10 -20 -30 -40 -50 
offset from HW2 (E-W) ["] 




50 40 30 20 10 -10 -20 -30 -40 -50 
offset from HW2 (E-W) ["] 



w 



4 * 



3 S 



-9.5 



-10 



■10.5 



-11 § 



■11.5 ' 



-12 



-12.5 



5 I 



Fig. 3. From top to bottom (a-c) the integrated flux, the veloc- 
ity field, and the line width, of the methanol 7_i — > 6_i E-type 
transition, the strongest unblended line in the observed spectra. 
In all maps the cross marks the "Centre" position, the diamond 
"NE Outflow" and the circle indicates the position where the 
"Envelope" spectra are evaluated. 



a secondary peak in the line width map. Throughout the paper 
we will refer to this position as the "NE outflow". The quiescent 
methanol emission that is not associated with the outflows and 
does not show any peak in the integrated line flux map we will 
refer to as the "Envelope". 



4. Analysis 

4.1. Rotation diagrams 

From the above spectral cuts in 3 positions and the brightness 
distribution of one of the lines, it is clear that the observations 
trace a complex region in which the dynamic and excitation char- 
acteristics change, at least on the scales represented by the beam 
size of the telescope. To make a representation of the large scale 
distribution and excitation of methanol, we have created rotation 
diagrams (Boltzmann plots) at every pixel in order to be able to 
study the characteristics of the excitation distribution. Inherent 
to the method are a number of assumptions, namely that the gas 
at each velocity and spatial resolution element can be described 
by a single excitation temperature (T lot = r ex ), that all the lines 
are optically thin, i.e. r « 1, and that the size of the emitting re- 
gion is the same for all lines. Although not all these assumptions 
may hold for all positions, the strength of the method is that it 
shows general trends on the relevant scales of our observations. 
In Sec. l4.4l we will test the validity of these assumptions. 

For each pixel at which we measure at least three methanol 
transitions with S/N>5, we plot the upper energy level E u [K] 
of the transitions on the x-axis versus the logarithm of the col- 
umn density in the upper energy state divided by the statisti- 
cal weight of that upper level log(N u /g u ) on the y-axis. The 
weighted column den sity (N u /g„) can be calculated with Eq. 14.11 
dHelmich et al.l [1994T) and the appropriate coefficients in Table 
Q] In the equation Q(T mt ) is the partition function, jj the perma- 
nent dipole mome nt [Debye], S g the line strength value from 
dBlake et alJI 19871) . and Tmb the main beam temperature is the 
antenna temperature scaled by the main beam efficiency. 



N M 
Q(T mt ) ( 



-EJT rc 



1.67 • 10 1 

VfXpg 



I 



T M ndv 



(1) 



By fitting a straight line through these data points we deter- 
mine the rotation temperature T ml [K] and column density Nm 
[cm -2 ] of the methanol gas. The rotation temperature is the neg- 
ative inverse of the slope of the fitted line and the total column 
density is where the extrapolated line crosses the y-axis. 

Fig. [4] shows the resulting rotation diagrams of the 
"Envelope" (the lower temperature quiescent gas), the "NE out- 
flow", and of the first velocity component at the "Centre" posi- 
tion of HW2. For the two positions in the outflow and envelope 
a single gas component provides a good fit to the data points 
with a rotation temperature of 30 K to 35 K and a column den- 
sity of 10 16 crrT 2 . However, in the centre region things are more 
complex. The velocity component at -10.5 km s~' can be fitted 
by a straight line, yielding a rotation temperature of 224 K and 
a column density of 2xl0 15 cirT 2 . For the -5.2kms _1 velocity 
component, the points lie nearly horizontal in the diagram and 
we would only able to put a lower limit of >300 K on the rota- 
tion temperature with the simple fit. 

4.2. Population diagram modeling 

Clearly the basic assumptions going into the rotation diagram 
analysis are violated for the "Centre" position and a successful 
explanation must accommodate both an apparent negative rota- 
tion temperature and the relative weakness of the low K lines for 
the red-shifted component. Non-LTE conditions, optical depth 
and beam-dilution effects can all play a role here. As a first 
step, to investigate the effect of line optical depth on the above 
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Fig. 4. Rotation diagrams for three specific components and a population diagram for the fourth component identified in Sec. I3.2I 
Triangles indicate upper limits and error bars the lcr standard deviation in the flux determination. The diagrams displayed are: a: 
rotation diagram for the main, -10.5 km s~' velocity component at the "Centre" position (+3", +3")- b: population diagram for the 
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includes optical depth effects and beam dilution (see Sec. 14. 21 , c: rotation diagram for the "NE Outflow" (+27", +9"), d: rotation 
diagram for "Envelope" (+21",-3"). 



results , we have adopted the method of Goldsmith & Langer 
(1999) which includes the source size (8 S ) as a free parameter 
and allows for optical depth effects. We have performed a x 1 
analysis for 10 < r rot < 500 K, 10 14 < Nm < 10 18 cnr 2 , and 
1" < 9 S < 14", i.e. a beam dilution factor between 200 and 1. 
For the red-shifted component, we find the best fit to be for a ro- 
tation temperature of 125 K, a column density of 4 x 10 17 cm -2 , 
and 9 S = 1", much smaller than the beam. The merit of this 
model can be seen in panel Fig.[4j3. In this specific case the max- 
imum optical depth encountered was r = 6.8. As this seems to 
be the component with the highest column density it is not sur- 
prising that 13 CH30H is detected in this component. Obviously, 
the highly excited K lines lead the model towards a high col- 
umn density, but at the same time, the only way for the model to 
accommodate the lower K lines by increasing the optical depth 
is to decrease the source size. It should however be noted that 
this analysis assumes a single excitation temperature. We will 
discuss non-LTE modeling in Sec. 14.41 

4.3. Spatial distribution of the excitation 

Maps of the derived rotation temperatures and column densities 
are presented in Figs. [5] and [6] In the maps we have excluded the 



results of the -5.2 km s velocity component seen at the position 
of HW2. The methanol rotation temperature ranges from ~20 K 
to ~100 K over most of the map, but shows a pronounced peak 
of >200 K close to the HW2 position. We find column densities 
between 10 15 cm -2 and 3 x 10 16 cm -2 . Overall, the methanol 
column density distribution is quite smooth, showing a mini- 
mum towards the centre position. In Sec. 14.41 we will investi- 
gate whether this minimum is an artifact of the rotation diagram 
method. Additionally, some care must be taken when analysing 
the results in the centre region. Only at the position of HW2 are 
we able to separate the two velocity components of the gas. In 
the neighbouring pixels the integrated line fluxes are therefore 
overestimated due to blending. 

4.4. Non-LTE analysis 
4.4.1. Method 

In order to check the validity of the assumptions of the rota- 
tion diagram analysis, we have run a set of non-LTE models. 
By comparing the resulting synthetic spectra with the observed 
spectra, these models also allow us to investigate the physical 
conditions (i.e. temperature and density) at distinct positions in 
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Fig. 5. Map of derived rotation temperatures. Symbols are as in 
FigE] 
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Fig. 6. Map of derived column densities. Symbols are as in Fig[3] 

the methanol emitting region. We have used the RADEXB pack- 
age, that uses statistical equilibrium for the population calcula- 
tions a nd an escape prob ability method to calculate the optical 
depth ( Ivan der Tak et al.ll2.007h . The main assumptions are that 
the medium is uniform and that the optical depth is not too high 
(< 100). The most accurate c ollisional rates are available for the 
first 100 levels of methanol dPottage et al.ll2004t ISchoier et ail 
12001 and we therefore restrict this part of the analysis to the 24 
relevant lines from Table Q] 

Initial runs of the non-LTE model indicated that the methanol 
excitation is more strongly governed by the density than the 
temperature. Using a uniform sphere geometry, we considered 
kinetic temperatures between 30 and 300 K and H2 densities 
between 10 4 and 10 7 cirT 3 . Even at the lowest temperature 
all the lines, except the K=+5 and K=±6, are excited in the 
high-density case. On the other hand, at low-density and high- 
temperature, only the lowest three K levels (K=0,+ 1,+2) are 
populated. A combination of high temperature and high density 
is required to excite all lines, maki ng this methanol ba nd useful 
as a density and temperature tracer dLeurini et al.ll2004l) . 

Starting with model values comparable to the results from 
Sec. 14. II it is possible to model the observed line strengths of the 
large-scale emission with the simple assumption that the emit- 

2 http://www.strw.leidenuniv.nl/~moldata/radex.html 



ting region is filling the beam. For these parameters all the lines 
are optically thin (r j 0.3), validating the assumption made in the 
rotation diagram analysis for the large scale emission. 

In the above calculations the background radiation field 
was represented by the cosmic microwave background with 
7bg=2.73 K. Models with background radiation temperatures in 
the range 3 - 300 K were run, but this has no clear impact on the 
lines considered here. Mostly the contrast of the lines with re- 
spect to the background was reduced without changing their rel- 
ative intensities. This result can be understood from the fact that 
our models consider only the to rsional- vibrational ground state 
of methanol dLeurini et al.l2007t) . These models are inconclusive 
about the influence of the radiation field, which are expected to 
be important in the maser region. 

4.4.2. Non-LTE results 

Next, we use the non-LTE calculations to fit the spectra for 
the four regions identified in Sec. 13.21 (both velocity compo- 
nents at "Centre", the "NE outflow", and the quiescent gas in 
the "Envelope"). In a search for minima of the reduced x 1 we 
have used a model grid that includes column densities from 
10 13 to 10 19 cm~ 2 at every half dec, densities from 10 4 up to 
10 9 cm 4 , kinetic temperatures between 30 - 300K, and radia- 
tion temperatures in the range 3 -400K. Based on our findings in 
Sec. l4.ll we have adopted a beam dilution factor of 50 for the red- 
shifted emission in the "Centre", corresponding to the 2" extent 
of the region where the masers appear. It is found that the mod- 
els are able to reproduce the lines at all four positions to within 
5-10% accuracy. In Fig. |7]we plot the observed integrated flux 
and the synthetic integrated flux calculated with RADEX versus 
7k — >6k number. Also indicated is our 3<x detection limit. 

Like before, this analysis is complicated somewhat by the 
blending of lines, both by methanol and in the case of the K=+l 
E by blending with SO2. In the case of the K=+l E line we 
have not included it in the analysis of the centre position com- 
ponents as the SO2 contributes significantly to it. For the "NE 
outflow" and for the "Envelope" the SO2 does not appear to 
contribute significantly, as the other SO2 line is not detected, 
and we have included it in the analysis. At the position of HW2 
the methanol gas is highly excited and the blending of methanol 
lines becomes more severe and confusion prevents one from de- 
tecting the higher (K> 5) lines. Also, at this position the K=+2 
lines show a large (~ 50%) deviation from the RADEX model. 
The reason for this discrepancy is yet unclear, but could be an- 
other indication that we are not complete in treating radiative 
processes. 

The results of our non-LTE analysis are summarised in Table 
[3] The coolest and most diffuse gas is found in the envelope 
and outflow, where one should take into account that the gran- 
ularity of the model grid probably prevents the distinction of 
small temperature and density differences. In comparison, the 
gas associated with the rotation temperature peak is much denser 
and also warm. We find the highest density and temperature 
in the "Centre" HW2 position where the maser emission oc- 
curs. However, given the large beam dilution and the limited set 
of available levels for collisional and radiative excitation, these 
numbers should be treated cautiously. 

The non-LTE analysis confirms that there are no large col- 
umn density gradients across the source. But although gener- 
ally the values are in agreement with our rotation diagram anal- 
ysis, they differ in details. In particular, the non-LTE analysis 
finds local maxima of column density in the "Centre" compo- 
nents. Although some shortcomings are demonstrated for the in- 
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Fig. 7. Plots of the integrated intensity for different CH3OH 7k — > 6k transitions and the best x 2 model values. In the panels 
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ner components, the rotation diagram maps are useful for out- 
lining the qualitative large-scale properties and constraining the 
distribution of the excitation; running pixel based non-LTE mod- 
els is beyond the scope of this paper. 

4.5. H 2 column density 

In order to estimate the methanol abundance we require a mea- 
sure of the H2 column density to combine with the estimates of 
the met hanol column density. There fore we used the SCUBA 
map of iBottinelli & WUUamsl d2004l) and derived the H 2 col- 
umn distribution using the formula of iHenning et al.l (l2000l) (eq. 



2), where S v is the flux in Jy beam -1 , Kd{A) the dust opacity at 
the observed wavelength, fi m /, the solid angle of the main beam, 
B V (T) the Planck function (Black body), mg the mass of the hy- 
drogen atom, and R is the gas to dust rati o. In the calculation w e 
have used: /c d (875//m) = 1.4 andfi=150 dHenning et alJl2000h . 



K d (A)Ll mb By(T)2m H 

The calculations were run for two cases. In the first case, we 
assume a constant dust temperature of T&= 30 K to derive the H2 
column density and the methanol abundance, Xm- In the second 
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Table 3. Results of the non-LTE models. 



Position 




AS* 


log(N M ) 


Tkin 


log(n H2 ) 


Trad 








-2 

cm - 


K 


cnT 3 


K 


Centre, main 


+3 


+3 


16.0 


90 


4.0 


30 


Centre, red-shifted b 


+3 


+3 


17.0 


130 


7.0 


100 


Outflow (NE) 


+27 


+9 


15.0 


90 


5.0 


30 


Envelope 


+21 


-3 


15.0 


90 


5.0 


30 



a Offset from centre of map RA 22' ! 56 m 17.88', DEC +62°01'49.2". 
b This component was modeled with a beam dilution factor of 50. 
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Fig. 8. H2 column densities derived from the SCUBA dust con- 
tinuum for a constant dust temperature of 30 K. 
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Fig. 9. Methanol abundance derived from the methanol column 
densities obtained in the rotation diagram analysis and the H2 
column densities obtained for the constant dust temperature 
30 K. 



case we have assumed the dust temperature to be equal to the 
rotation temperature that we derived from the rotation diagram 
analysis in Sec. I4.ll and used that as input to the black body 
function. 

4.6. Methanol abundance 

The results of the constant temperature hydrogen column density 
and methanol abundance calculations are presented in Figs. [8]|9] 



The res ulting hydrogen column dens ity Fig. [8] is the SCUBA 
map (by Botti nelli & Williams! (12004)) scaled by a constant fac- 
tor to convert it from Jy beam 1 to N(H2) in cm -2 . The re- 
sult is a very smooth map with H2 column densities between 
~ 10 21 cm~ 2 and ~ 10 24 cm 2 . The general morphology of the 
dust emission is similar to that of the methanol gas, though no- 
tably the dust emission does not peak at the position of HW2 but 
further to the SE at position (-9",-4"). 

Similarly, the methanol abundance map derived with a con- 
stant dust temperature, Fig. [9] shows the same structure as 
the column density map, since there are no large variations 
in the SCUBA dust map. We derive a methanol abundance of 
Xm ~10~ 85 for the central region. As noted before, the distri- 
bution is relatively flat, and probably not very accurate for the 
centre position. 

In the second case, when we used the methanol rotation tem- 
perature as a measure of the dust temperature, we find the result- 
ing H2 column density to be anti-correlated to the rotation tem- 
perature, as the high rotation temperatures derived for the centre 
region then result in a lower H2 column. In this way we would 
derive a H2 column of ~ 10 22 cm" 2 (much lowe r than inferred 
earlier ~ 10 24 cm' 2 ; iMartfn-Pintado et al.ll2005b . Although we 
have already encountered the limited validity of the rotation di- 
agram for the "Centre" region in Sec. 14.4.21 it is clear that there 
must be a mismatch between the gas and dust temperature on the 
scales of the JCMT maps. Both the rotation diagram and non- 
LTE calculations indicate higher temperatures towards the cen- 
tre, for which there is little evidence in the dust images, which 
even peak at an offset of ~10" from the rotation temperature 
maximum. Besides the possibility that the rotation temperature 
is not a good estimate of the kinetic temperature of the gas, it 
seems plausible that the kinetic temperature of the gas is not 
equal to the dust temperature at these locations. Therefore, we 
adopt a constant dust temperature of 30 K to estimate the H2 col- 
umn density and CH3OH abundance, realizing that this method 
will overestimate NiHj) at the centre of the map. 

5. Discussion 

5.1. Outflow morphology 

We have presented maps of the large scale distribution of the 
thermal methanol associated with Cep A East (Figs. [3] and 
[6). The emission is constrained to a region of size 0.29 x 
0. 15 pc (85" x 43"), roughly centred on the location of the radio 
continuum source HW2, the most massive YSO in the region. 
Moreover, the velocity field shows that the methanol gas is en- 
trained in a large scale bipolar outflow. The position angle of this 
outflow appears consistent with the HW2 geometry and earlier 
results (iGomez et al.ll 19991) . 

While the largest flux of the K=-l E line Fig. [3^ is associated 
with the receding part of the outflow (towards the SW), coinci- 
dent with the SCUBA peak, the most energetic gas appears at 
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the opposite side of HW2. This is evident already in the velocity 
width distribution of the brightest line, that has a line width of 
>4 kms" 1 , to the NE of HW2, Fig. [5J;. It becomes more clear 
when the temperatures are derived from rotation diagram analy- 
sis, the rotation temperature peaks to the NE of HW2 at ~200 K 
(Fig.0. 

Several authors dPatel et al.l2005t Uimenez-Serra et al.ll2007t 
iTorrelles et al.ll2007l) have inferred obscuring gas in a flattened 
circumstellar structure of ~1" (700 AU) surrounding HW2. 
The gas shows signs of rotation and is perpendicular to the 
high velocity outflow ( ~500 kms" 1 ) seen in the radio contin- 
uum (ICuriel et al.ll2 006). The molecular outflow that we observe 
in methanol is blue-shifted to the NE, which implies that this 
side of the outflow is p ointing towards the observer, Fig. [3j>. 
iTorstensson et al.l (1201 ll) propose a model in which the methanol 
maser emission arises in a ring-structure in the equatorial plane 
of HW2. In contrast to the circumstellar molecular material, 
the masers do not show signs of rotation; rather the observa- 
tions seem to indicate an infall signature. It is argued that the 
masers occur in the shock interface between the circumstel- 
lar disk and inflow reg u lated by a large scale magnetic field 
dVlem mings et al. 2010). ITorstensson et al.l d201 ll) derive an in- 
clination of 67.5°and a position angle of 9.3° for the ring. The 
masering gas occurs on size sc ales of ~2" and is very likely 
clumpy on smaller scales ~0.1"(Minie r et al. 1 12001 . At the res- 
olution of the observations presented in this paper the masering 
gas can clearly not be resolved, nevertheless the geometry that 
is observed on smaller scales suggest that to the NE of HW2 we 
have a line of sight into the outflow probing highly excited gas 
close to the protostar, which explains why we see such a high ro- 
tation temperature in this region. Evidence for s uch wide-angle 
molec ular outflow has recently been presented bv lTorrelles et al.l 
(120101) . 

5.2. Methanol distribution 

The highest methanol temperatures appear to occur close to the 
HW2 central source, but we have demonstrated in Sec. I4.6l that 
the dust emission is not showing the same distribution. It seems 
plausible that the highest methanol rotation temperatures occur 
in a region where the dust is at least partly destroyed, possibly 
at the onset of the outflow where the evaporation from the grains 
occurs that produces the methanol gas in the first place. The 
considerations on the dust distribution force us to adopt a con- 
stant dust temperature in our abundance analysis. On small radii 
this may have two effects: firstly, the gas to dust ratio may be 
higher due to grain destruction, which increases N(H2) and low- 
ers the abundance. Secondly, the dust temperature may be higher 
than 30 K, which decreases NiHi) and increases the abundance. 
Increasing the gas-to-dust ratio from 150 to 200/250 decreases 
the abundance with a factor of 0.75/0.60 respectively. Increasing 
the dust temperature to 60/100/300 K increases the resulting 
abundance by a factor of 2.3/4.1/13.0 respectively. High optical 
depth of the dust will lower the derived abundance. These argu- 
ments show that at the peak rotation temperature, the methanol 
abundance is quite sensitive to increased dust temperatures, but 
is probably correct to within an order of magnitude, without tak- 
ing beam-dilution effects into account. 

The derived abundances of ~ 10~ 8 ' 5 are higher than ex- 
pecte d from gas-phase chemistry, but lower than seen in the solid 
state dGibb et alJl2004l) . High optical depth of the lines is ruled 
out by the excitation calculations, except in the central region. 
However, it is possible that there the local methanol abundance 
is still higher. For example, the methanol maser occurs on much 



smaller scales than those we are probing in these JCMT observa- 
tions. The total extent of the CH3OH maser is ~<2". Moreover, 
we found that the same small scale enhancement can explain 
the peculiar excitation of the "Centre" region. This strengthens 
the idea that the highly excited methanol gas emission is com- 
ing from the same region. Assuming that the dust and the cor- 
responding H2 is smoothly distributed, this component can have 
higher column density because of the optical depth and beam- 
dilution, resulting in methanol abundance estimate of ~ 10~ 7 lo- 
cally. This valu e is similar to previous valu es by other authors for 
maser regions dMenten et al] 1 19861 1 19881) and still higher than 
expected from gas-phase chemistry and lower than seen in the 
solid state. Possibly the methanol evaporated from grain man- 
tles some time ago (~ 1000 yr) and is now destroyed by chemical 
reactions and/or photo-dissociated. 

As noted earlier, all the methanol gas seem to be en- 
trained in an outflow, spanning 0.29 pc. The outflow is ori- 
ented clo se to the plane of the sky, assuming an inclination 
of 67.5° (Tors tensson et al.l I2.Q1 ll) we find an outflow velocity 
of ~5.8 kms" 1 . The result implies a dynamical timescale of 
~ 2.4 x 10 4 yrs. The dynamical timescale agrees well with the 
chemical timescale (a few 10 4 years) fo r which such a methano l 
enhancement is believed to be possible (Ivan der Tak et al . 2000). 
This supports the argument of a common driving source for 
all the methanol gas in the region. Alternatively, the methanol 
may be released in the wake/ outflow cavity of the jets described 
bv ICunnin gham et al. (2009): the uv-radiation produced in the 
shock interface may shine back into the cavity and photo-desorb 
the methanol, causing a local enhancement in the methanol col- 
umn density. 

5.3. Physical conditions 

The rotation diagram analysis shows rather warm gas entrained 
in a large scale outflow. Although the non-LTE analysis demon- 
strates that the methanol band used in this study is quite sensitive 
to collisional excitation, even in the low temperature regime, it 
confirms qualitatively the large scale presence of warm methanol 
gas. Over large areas of the source the total methanol column 
density appears to be fairly constant at ~ 10 15 cm" 2 , assuming a 
uniform distribution with unity beam-filling factor. 

The fidelity of the rotation diagrams is breaking down par- 
ticularly in the centre area, where higher excitation lines ap- 
pear. Especially the red-shifted component requires non-LTE 
excitation and/or beam-dilution effects to be explained. At 
these locations a complete analysis requires high resolution 
data and non-LTE modeling including excitation by infrared 
radiation. That more methanol levels should be taken into ac- 
count is demonstrated by the detection of the torsionally ex- 
cited - 1 A v t =l line at 33 7.969 GHz in the red-shifted com- 
ponent dLeurini et al 120071) . Radiative pumpi ng is also expected 
from models of methano l maser excitation (|Cragg et al.l 12005 ; 
ISobolev & Deguchill 1 994t ISobolev et all 199% 

5.4. Origin of maser emission 

Although a larg er area around Cep A has been searched for 
maser emission dTorstensson et al.ll201 ll) . only one place in this 
source, so abundant in methanol, does show maser emission. We 
find the maser position to contain warmer gas than the other 
positions, probably the highest column density, and certainly 
the density appears to have a maximum here at ~ 10 7 cm" 3 . 
Moreover, this position does stand out as the position where a 
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second kinematic component at -5.2kms _1 is found. Special 
properties for this component are the detection of the torsion- 
ally excited line, and the requirement to invoke a large beam 
dilution to explain the ratios of the higher excited lines. This 
dilution factor matches quite well with the size of the area of 
< 2" over which mas ers are detected (Sugivama et al. 2008; 
Torstensson et al. 201 1). Combining these facts, we identify this 
component with the origin of the maser action. We note that the 
densities we derive for this component are still below 10 9 cm -3 , 
where the m aser lines would b e quenched according to model 
calculations (Cragg et al. 2005). The models have methanol tem- 
perature ranges and methanol column densities that are consis- 
tent with our estimates, but it is hard to test the dust temperature, 
which must be >100 K for the maser excitation. We note that 
many of the values estimated here are inferred indirectly and 
the relevant processes may occur on even smaller scales than re- 
flected by the beam-dilution factor. 

The central abundance of Xm ~ 10 -7 is not in excess to that 
found in low-mass star-forming regions, bu t the methanol col - 
umn density N M ~ 10 17 cirT 2 clearly is (Ma ret et al.1 12005). 
This could point to the condition that the high-density, high- 
temperature region in low-mass star-forming regions is too small 
to produce a path length long enough, with suffi cient methanol 
colum n density, to produce maser emission. Ij0rgensen et al.l 
(2002) show the high-density, high-temperature region for a low- 
mass protostar to be ~ 10-30 AU. In contras t, for a high-mas s 
protostar, the region extends over -1000 AU (iDotv et al.ll2002l) . 
possibly facil i tating the necessary amplification. Alternatively, 
Pandi an et all |2008) suggest that in low-mass stars the critical 
infrared pumping is only available in the inner region, where the 
density is so high that the maser will be quenched. Both effects 
favour maser emission in the environments of high mass proto- 
stars and explain that th e masers occur on rathe r large distances 
from the central object (Torstensson et al. 20TlT). 

To further constrain the physical conditions at the position 
of the methanol maser and discern radiation and density effects, 
high-J and v, = 1 line observations are required, for example 
with the CHAMP+ instrument on the APEX telescope. To probe 
the excitation of methanol gas on small scales (comparable to 
that of the maser emitting region) high resolution interferome- 
try observations are needed, with instruments such as the SMA 
and the IRAM interferometers, and in the longer term, for lower 
declination sources, ALMA. The interpretation of these observa- 
tions requires models that include the excitation of CH3OH by 
infrared radiation. 
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